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Abstract 



Classical defects (monopoles, vortices, etc.) are a characteristic conse- 
quence of many phase transitions of quantum fields. We show a model in 
which the onset of classical probability distributions, for the long- wavelength 
modes at early times, allows the identification of line-zeroes of the field with 
vortex separation. We obtain a refined version of Kibble's causal results for 
defect separation, but from a completely different approach. It is appar- 
ent that vortices are not created from thermal fluctuations in the Ginzburg 
regime. 
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Because phase transitions take place in a finite time, causality requires correlation lengths 
to remain finite. As a result, scalar order parameter fields 0(x) become frustrated, and 
topological defects arise so as to reconcile field phases between different correlation volumes 
mj^. For example, the breaking of the GUT symmetry in the early universe gives rise to 
monopoles, the most familiar defect. However, they are but one of many possibilities, that 
include cosmic strings (vortices), which may be the source of high energy cosmic rays, as 
well as contributing to structure formation in the early universe. 

For the weak-coupling theories of the early universe, defects have non-perturbatively 
large masses, comparable to the temperature scale ksT^ at which the transition takes place. 
Thus monopoles produced in the early universe at a GUT transition might be expected to 
have a mass m ~ lO^^GeV. For this reason alone, they and other defects are manifestly 
classical, even though the phase transition that produced them is intrinsically quantum 
mechanical. In this letter we show how the creation of classical defects comes about for 
simple temperature quenches with a finite quench rate through the critical temperature T^. 
[The case of an instantaneous quench has been considered elsewhere 0]. Finite quench rates 
have been studied by Kibble [Q, using causality bounds, from which he has derived simple 
scaling laws that incorporate mean-field dimensional analysis. 

Although the equilibrium correlation length ^eq(^) diverges at the transition, the true 
correlation length ^{t) does not. Thus, causality imposes a maximum rate at which the 
correlation length can grow and hence a maximum correlation length, C,, at the onset of 
the transition. At the same time, causality imposes a horizon outside which the fields are 
uncorrelated. If we simply relate the correlation length ^ and the separation length between 
vortices ^def by ^ = O(^def), then the density of defects is bounded, and calculable, at their 
moment of formation. 

If, in the vicinity of T^, dT/dt = —Tc/tq, then dimensional analysis suggests that the 
earliest possible time at which defects could be formed after the onset of the transition 
is t = O(tK), where = (tq//^^)^^^- This is deduced from .^eq(^ = ^k) = —1, where 
^eq(^) = |^~^(^)| = Ai~"^-^TQ/t, (/i is the particle mass). We shall use this as a benchmark in 
our calculations. Therefore, from causality arguments, we derive the bound |^ 

= at = tK) = u ~ M-^(rQ^)^/^ (1) 

leading to the estimate ndef = 0{^j^'^) for the defect density at the time of their production. 

But there is an alternative scenario for counting defects after the transition is completed. 
Simple defects have false ground-state or vacuum at their cores where the field vanishes. 
Under suitable conditions, the separation length for defects is more sensibly derived by 
counting zeroes of the field as ^def = O(^zero), where ^zero measures the separation of field 
zeroes. However, in principle, ^zero and ^def are different correlation lengths, since not all 
zeroes are candidate defects (because zeroes occur on all scales). One needs to count zeroes 
of an appropriately coarse-grained field, in which structure on a scale smaller than an specific 
classical vortex size ~ f^~^ = ^5 is not present. At a given time t* the separation of zeroes 
is ^zeroit*)- Then, in order that line-zeroes can be identified with classical defects vortex 
cores, and ^^emit*) with ^def we need satisfy, essentially, two important conditions 1^: I) the 
separation between zeroes must be insensitive to the cutoff scale A~^ ~ and II) the 
energy in field gradients should be commensurate with the energy in classical defects with 
the same density as that of line-zeroes. 
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We will show how decoherence can help to satisfy these two conditions. For weak cou- 
plings (massive defects) our main prediction agrees qualitatively with the causal (dimen- 
sional) prediction for counting defects, but here from a completely different point of view, 
in which classicalisation of long-wavelength (unstable) modes of the field warrants the iden- 
tification of line- zeroes with defects at the time in which the transition has been completed. 

There are further complications according as the symmetries are global or gauged 0, 
but the simplest of all are global, and it is these that we shall consider here. We have global 
vortices in mind, the simplest sensible defects, in which case the order parameter field $(x) 
is a complex scalar $ = (0i -|- i02)/v^, with action 



With /i^ > 0, the 0(2) symmetry is broken at the scale |$| = rj^rf = /i^/2A. [Generalisation 
to global monopoles, or even domain walls, is straightforward.] 
Global vortices are line defects in the field, solutions to 



around which the field phase 9 = he*^) changes by 27r. Considered as tubes of 'false' 
vacuum, with cold thickness 0(/i~^), they have a line of field zeroes at their cores and, 
once the transition is complete, energy per unit length a = Oi^jj^/X) (up to multiplicative 
logarithmic terms ln(yU,^def))- 

There is no unambiguous definition of the onset of classical behaviour, least of all for the 
production of classical defects. Our approach stresses the role of classical stochastic field 
equations. They arise in two forms, in Lagrange's form as stationary phase approximations 
to the evolution kernel of the density matrix, or in Hamilton's form as the dominant path in 
phase space of the evolution kernel of the Wigner functional. Initially we shall concentrate on 
the former. In general, if the system is closed, quantum interference between different field 
configurations forbids us from identifying the dynamical solutions to (^ for the evolving 
system as describing real defects. However, if the system is open, the environment with 
which it interacts can eliminate the interference, in which case the system can be said to 
have decohered. 

In practice, the systems of interest to us are open. This is true, firstly, in the trivial 
sense that, in the absence of superselection rules, our order parameter fields $(x) interact 
with everything else in the universe (termed x(x))) but this environment is traced over when 
looking at the transition. Further, for the continuous transitions of (|^) that will interest 
us here, field ordering takes place through the exponential growth of the unstable {k < n) 
long-wavelength modes $<(x) of $(x) . The stable {k > /i) short-wavelength modes $>(x) 
do not become classical and behave as a further part of the environment in which the long- 
wavelength modes, which constitute our 'system', decohere. As a consequence, the interiors 
of vortices are not classical. Since it is the field profile in the interior of the classical vortex 
that carries the non-perturbatively large energy or tension, this might be thought to be 
damaging, but we will see that it is not the case. In particular, topological charge can be 
classical. As a result, non-classical cores are irrelevant as far as counting defects is concerned. 
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Our classical equations are, most simply, stationary phase approximations to the evolu- 
tion of the reduced density matrix Pr[0o 4><yt] = (0<|Pr(^)|0<), obtained by tracing out the 
environments x ^i-nd 0> (we are using 0> to denote the real Cartesian doublet) . In a field 
basis (driven by our need to describe defects through field equations) Pr[0<, 0<, t] evolves in 
time by means of the coarse-grained effective action (CGEA) y4[0^,0^], as we have shown 
in Refs. (and in a more detailed version in 0). The CGEA takes into account the effect 
of the environment on the dynamics of the system. We regard the imaginary part of A as 
coming from a classical noise source ^(x) (or several such sources), say. If its probability 
distribution is p[^], then the CGEA can be rewritten as [|,|T3] 

40+0-] = -^In/ Pep[e]exp|25eff[0+,0<,e]}, (4) 

where 

5eff[0t'0<,e] = ReA[0+,0-] - I d^A{x)ax)], (5) 

and A(x) depends on the way P, pT| , p!^ in which the environment couples to the 0< field. 
For the simplest case of biquadratic couplings (e.g. 0<x^ or the inevitable 0<0>), we find 
A = (|<l>+|2- |$-|2)/2. 

The stationary phase equation, obtained by taking the functional variation 



SScs[(t>t'4><,^] 



5ReA[0+,0-,e] 



0, (6) 



is a generalised Langevin equation, with classical noise ^ and dissipation, to satisfy the 
fiuctuation-dissipation theorem (see Ref. for a complete deduction of Eq.(^). If we look 
for vortex solutions $< to (^, modelled on the vortex solutions to (|^), we shall find ourselves 
in some difficulty at early times, in that the noise will cause fiuctuations (e.g. small vortex 
loops) that will make the identification of a single vortex difficult. 

However, the existence of classical defects depends on more than clean vortex solutions 
to Eq.(|^). The absence of quantum interference between the stationary phase solutions 0< 
to the classical stochastic equations and their neighbours is manifest through the increasing 
diagonalisation of Pr[0o0<;^]- This leads to the crucial notion of a decoherence time to, 
after which Pr[0o '^<' ^] (°^' viol's exactly, its real part) is effectively diagonal. The equation 
(^) is the classical stochastic equation that we are looking for but, as it stands, is only 
guaranteed to describe classical defects after decoherence. 

A digression is necessary. There is not a universal decoherence time t^. toik) depends on 
wavelength. Long- wavelength field modes decohere first, shorter ones later, and those with 
k > n never. To make the problem of diagonalisation tractable we have restricted ourselves 
to field configurations $+(/co) peaked around a particular wavenumber [0,0^ • The order 
parameter itself corresponds to /cq = but, more generally, we are motivated by the way in 
which the power in the field fiuctuations is increasingly peaked around long-wavelengths in 
the early stage of the transition P, |T3tp!4| . That is, if the power spectrum P{k,t) of the field 



fiuctuations is defined by 

G<(r,t) = ($<(x)<|.*<(0)), = / ^,P{k,t)e^'^- (7) 
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then, provided the quench is fast enough, k'^P{k,t) rapidly develops a 'Bragg' peak at 
k"^ = kp(t) < iJ?. A sufficient condition for this to happen is that /irg ^ rj/fi [|14|. Tighter, 



but less transparent bounds can be given |]I4[. We restrict ourselves to quenches permitting 
such a dominant momentum peak and identify fco with it. The relevance of this is that it 
is the peak in k'^P{k,t) that sets the scale for the separation of line-zeroes in the field (j). 
Once kp{t) <^ it does not matter where the boundary A = 0{fi) between system and 
environment is set. Since defects can be characterised by the line- zeroes at their cores |]T3|, it 



is the dominant wavenumber feg that, after the onset of classical behaviour, sets the density 
of the vortices that we wish to determine. This peak at /cq also warrants the fulfilment of 
condition (I). 

This peaking stops once the field has sampled the ground state values at |$| = 77. It 
does this by the spinodal time tgp, defined as that value of t (after the critical temperature 
Tc has been crossed) for which = 77^. For the critical temperature T^. traversed at a 

finite rate Tq^ (/itq 3> 1) we find P, P^JT^ 



ewHr-^) \-^. (8) 




In Eq.(|D the time is the same Kibble causal time introduced earlier (before Eq.(|I])). 
Because of the exponential nature of the relationships, in practice tsp = ©(tx); but larger 
than it. 

For t ^ tsp, perturbation theory can be used to calculate to. Decoherence at scale ko 
occurs when the non-diagonal elements of the reduced density matrix are much smaller than 
the diagonal ones. In we have evaluated tjj for the longest-wavelength modes, and found 
that t£,{kQ = 0) < tsp for all sufficiently fast tq for a realistic class of environments (for the 
instantaneous quench see Refs. We have considered an explicit environment consisting 

of a large number of weakly coupled x fields, with coupling constants gi satisfying Y,i df = 
O(A^), which permit a consistent loop expansion, with /i <C Tc <C ?], and we continue with 
this hierarchy here. The expressions are complicated, even in one-loop, and the reader is 
referred there for details. The main result for the decoherence time in a slow quench (and 
ko ^ 0) is 



3n 2 



4 ftD\ 2 1 r] 



Nevertheless, what interests us here is the shortest-wavelengths that will have decohered 
by time tgp. Provided ko is not too close to /i it is sufficient to replace nt^ by t£){ko)\J fj? — k^ 
in the calculations above to determine the bound on the decoherence time for the mode kQ. 
An immediate question is to ask what is the maximum value /Cmax for which tD(^max) = ^sp- 
After this replacement in Eqs.(^ and (^ we obtain 

72 / 2/3 ^ 



"max 



In 



1 yr,l/3 

Alternatively, there is a time t{kp) when the mode kp would become the dominant 
wavenumber. This 'quantum time' is determined by the instabilities of the long- wavelength 
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modes, which increasingly bunch about the wavenumber kp{t) ~ fi/t{kp) [Q. At time 



In . (10) 



Thus, the dominant modes at time tgp have aheady decohered by this time. 
Therefore, from the two previous equations we can see that 



^-2,*Kln(^J>l. (11) 

That is, the modes that determine defect separation are well decohered by time tgp- 
However, toik) oo as A; — » and the modes that characterise the classical profiles of 
individual defects (A; > jj) are not decohered. 

This makes constructing exemplary solutions to Eq.(^, and averaging them, a suspect 
activity. Rather, we examine the second strand of classicality, that of classical probabilities, 
as expressed through classical correlations. This is understood as the quantum system, in 
full (field) phase-space, mimicking the classical dynamics. The comparison is made through 
the reduced Wigner functional, defined as 



W\ 



As is well known, iy[0<,7r<;t] itself is not guaranteed to be positive until decoherence has 
occurred Jl^, although 



P[0<]i = j I?7r< W^[0<,7r<;t] = (0< |p.(t) |0<) > 0, (13) 

is always a true probability distribution for the field configuration 0< at time t. 

After decoherence, when Vr[0<,7r<;t] is positive, at least for long- wavelength modes, 
it can be identified with the Fokker-Planck probability distribution function P/^^[0<, 7r<], 
and P[(j)<^t can be equally identified as the Fokker-Planck probability P^^[0<]f. Since the 
Fokker-Planck equation is the obverse of the Langevin equation, at this time we can equally 
average over Langevin noise or over P^^[(j)^f 

However, as we observed earlier, the validity of the classical stochastic equations is, of 
itself, not sufficient to guarantee classical defects with energy profiles derived from (^). We 
do not have defects until the fiuctuations around them are negligible, and this is associated 
with the peaking of the power spectrum (and the validity of our approximation of adopting 
a dominant wavenumber) rather than decoherence or classical correlations directly |]5|,p!8|. 
To see this, suppose we can calculate P[(t><\t for all t. This permits us to calculate the equal 
time n-point correlation functions, in which we restore the subscripts a,b, ... = 1,2, 



'-^<ab..c 



Xi, .., Xn, t) = J P0<Pt[0<]0<a(xi)...(/)<c(Xn) (14) 



for all times t < tgp- As we know, equal-time correlators are all we need to calculate densities 
of field zeroes (line-zeroes, etc.). However, if all we are going to use is P[0<]t, the diagonal 
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matrix element of Pr(^), there is no real need to construct the Wigner functional. We can just 
do a calculation of P[0<]t from the start, along the lines of Ref. The weak interaction 



with the environment that is so successful in diagonalising the density matrix has only a 
small effect on the diagonal matrix elements P[0<]j for most of the time until tgp- 

The major characteristic of defects is their topological charge. We now see how localised 
topological charge precedes the appearance of vortices, and gives an estimate of their density 
when they do appear. Let us now consider the line-zero ensemble density Uzeroit) and 
the line- zero separation ^zero(^) for the long- wavelength mode fields. When the Gaussian 
approximation |II?|,pU[ is satisfied (as happens for t < tgp) they are determined completely 



by the short- distance behaviour of G^{r,t) as 

Since G<(r, t) has short- wavelength modes on the scale of a classical vortex removed, G'<(0, t) 
is finite. Initially, with P{k, t) large for k ^ n, G^{r, t) is very dependent on the value of the 
cut-off A. As a result, line- zeroes are extremely fractal, with a separation ^zero(^) proportional 
to the scale at which they are viewed, and are certainly not candidates for defects. Only 
once the peak a.t k = ko is firmly in the interval k < fi does ^zero(^) becomes insensitive to a 
cut-off 0(/i~^). This means that line-zeroes are straight at this scale, although they can be 
approximately random walks at much larger scales (condition (I)). 

Further, for this early period, when the self-consistent mean-field approximation is valid, 
the field energy {E)t of the system field 0< in a box of volume V becomes P,p!8[ 



(E), = V[(|V0<r), + A(r^2-(|0<r),n 

= V[27rn,Ut)G<{0,t) + \{r,^ - G'<(0,t))2] (16) 
= 27rL,ero(t)G<(0,t) + V[A(r/2 _ G^{Q,t)f]. (17) 

As before, x-field fluctuations are absorbed in the definition of yU^ = Xrf. Eq.(|l^) is obtained 
by using (|T5|) and L^croif) = Vn^ero is the total length of line-zeroes, on a scale in the 
box of volume V. 

We understand Eg. (p!7D , valid from time t = 0, when G(0, 0) = 0(/i^), until time t ^ tgp, 
when G(0, tsp) — rj^-, as follows. At early times most of the system field energy (proportional 
to V) is in fluctuations not associated with line zeroes, arising from the field potential. As 
time passes their energy density decreases as the system field approaches its post-transition 
value, becoming approximately zero. In addition there is a term, arising from the field gra- 
dients, proportional to the length L^ero of line-zeroes, whose energy per unit length increases 
from 0(/i^) to rf. At time tsp, when the fluctuation energy can be ignored, we have 

{E)t ~ Lzeroa, (18) 

essentially the energy required to produce a classical vortex tangle of length Lzero (up to 
0(1) factors from the logarithmic tails). Eq.(|18[) completes condition (II) in order to relate 
line-zeroes with separation length between defects. Although these line zeroes have the 
topological charge and the energy of classical vortices, they are not yet fully-fledged defects. 
However, by the end of the linear regime, when decoherence on the scale ^zero(^) has been 
effected, the final coupling of radial to angular modes that turns these proto-vortices into 
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vortices incurs no significant energy change, and n^evo of Eq.([T5D is a reliable guide for the 
initial vortex density. We stress that, as far as counting vortices is concerned, all that matters 
is how the power in the field fluctuations is distributed. The distance between defects is the 
relevant wavelength, and not the defect size, which shows no decoherence. 

This distance can be computed from the two point function of the field (see Eq.(|l5|)). 
Up to logarithms in tq and the other parameters of the theory it is given by p|, p!4| , p^ , 

~ Up) = W ~ Ai"^(T■Q/u)^^^ (19) 

in accordance with the predictions of Kibble 0, but for very different reasons. Kibble 
correctly argues that the long-distance adiabatic correlation length ^eq (that would diverge 
at T = T(.) freezes in at time tx before the transition with value ^, and this sets the scale 
for defect separation. However, in our picture defects do not form until time tgp after the 
transition with a separation given by that of line- zeroes, that is arbitrarily small at time ty^, 
and at the transition itself. That this is directly related to ^ is a consequence of dimensional 
analysis. This constraint was used to bound the production of monopoles and cosmic 
strings (vortices) in the early universe. [If we had broken an 0(3) symmetry (a = 
1,2,3) the corresponding defects would be global monopoles, giving qualitatively similar 
conclusions.] This is all for flat space-time. 

In summary, the mechanism for the production of classical vortices that we have pro- 
posed here has several parts. Firstly, the environment renders the long-wavelength modes, 
wavenumber k^, of the order-parameter field classical at early times tD{ko), by or before 
the transition is complete at time tgp- In particular, those on the scale of the separation of 
the line-zeroes that will characterise the classical defects will have decohered, even though 
the field modes on the scale of classical vortex thickness do not decohere. For all that, the 
field possesses classical correlations at early times by virtue of the quasi- Gaussian nature 
of the regime. For the longer-wavelength modes the classical behaviour of the field is ex- 
pressed through the classical Langevin stochastic equations that it satisfies. However, for 
line-zeroes to mature into vortex cores classical stochastic equations @ or, equivalently, 
classical Fokker-Planck equations are not enough. The field needs to have an energy profile 
commensurate with the vortex solutions to the ordinary classical Euler-Lagrange equations 
d). This requires that the field fluctuations are peaked around long-wavelengths, to avoid 
fluctuations causing wiggles in the cores and creating small cortex loops, a related condition 
satisfied in our models. The resultant density of line-zeroes can already be inferred in the 
linear regime, whose topological charges are well-defined even though close inspection of 
their interior structure is not permitted classically. 

The final result for their density, at the time tgp of their formation is, up to logarithms, 
that proposed by Kibble [Q from dimensional analysis of mean-field theory, although the 
reasoning is very different. Although logarithms introduce different scales in principle, in 
practice they have no qualitative effect. It is apparent that vortices are not created from 
thermal fluctuations in the Ginzburg regime, as suggested in the first instance by Kibble in 
an earlier paper 
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